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We discuss the classical and quantum mechanical evolution of systems described 
by a Hamiltonian that is a function of a solvable one, both classically and quan- 
tum mechanically. The case in which the solvable Hamiltonian corresponds to the 
harmonic oscillator is emphasized. We show that, in spite of the similarities at the 
classical level, the quantum evolution is very different. In particular, this difference 
q-i is important in constructing coherent states, which is impossible in most cases. The 

class of Hamiltonians we consider is interesting due to its pedagogical value and 
its applicability to some open research problems in quantum optics and quantum 
■ gravity. 

in 
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O ■ I. INTRODUCTION 



The goal of this article is to discuss the classical and quantum mechanics of systems 



whose Hamiltonian H = f(H ) is a function of the harmonic oscillator Hamiltonian H . 
The results can be easily generalized to other choices of H for which the classical and 
quantum equations of motion are exactly solvable. 

Once we solve the classical equations of motion for H , it is possible to study a system 
described by H = f(Ho). Although the solution is a straightforward exercise in classical 
mechanics, we will discuss it in detail because it is interesting to compare the solutions 
corresponding to both classical and quantum dynamics. Quantum mechanically the unitary 
evolution operator for H = f(H ) can also be constructed exactly once we know the evolution 
generated by Hq. A comparison of the dynamics given by Hq and /(-Ho) will allow us to 
analyze some distinctive features of the coherent states of the harmonic oscillator and discuss 
the difficulties that appear when we try to construct similar states for the dynamics generated 
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by f(Ho). This comparison will help us understand some aspects of the open problem of 
building appropriate semiclassical states for general Hamiltonians. 

Non-trivial systems whose evolution can be solved exactly both classically and quantum 
mechanically are rare. Usually, realistic systems are described by Hamiltonians of the form 
H = H + Hi, where H is a solvable Hamiltonian and Hi represents a perturbation. In 
most cases it is impossible to give the solutions to the equations defined H, so it is necessary 
to resort to approximation methods. The starting point of perturbation theory is the known 
dynamics generated by Hq. The simplest choice of Ho is the Hamiltonian of a free particle. 
However, if we are considering a system that has bound states, it is much better to consider 
a solvable H Q with bound states, such as the harmonic oscillator. 

In this paper we consider a different way to perturb a solvable Hamiltonian Hq by con- 
sidering a function of it. If this function is close to the identity, it will be possible to treat 
the system as a perturbation of Hq in the usual sense; otherwise it will provide different but 
still solvable dynamics. 

We point out that these kind of Hamiltonians appear in some physical applications, 
for example, in the context of quantum optics and classical and quantum gravity. For 
instance, the propagation of light in Kerr mediap], H| — media with a refractive index with 
a component that depends on the intensity of the propagating electric field — is described 
(for a single mode field given by the creation and annihilation operators a) and a and in the 
low loss approximation) by 

H = X (N 2 -N) = X :N^. = X a^a 2 , (1.1) 

where x is related to the susceptibility of the medium, and the Hamiltonian is a function of 
the number operator N = a)a. The symbol : : denotes normal ordering (creation operators 
to the left of the annihilation ones) and the operators a and a) satisfy the usual commutation 
relation [a, a'] = 1. 

Another situation where we find this kind of Hamiltonian is in general relativity. Einstein- 
Rosen [3| waves are cylindrically symmetric solutions to the Einstein equations in vacuum. 
The Hamiltonian that describes this system is [3, [5| 

If = 2(l-exp(-# Q /2)), (1.2) 

where Hq is a free (and easily solvable) Hamiltonian. 

The examples we have mentioned are field theories with Hamiltonians f{Ho) that are 
functions of free Hamiltonians Hq (i.e. quadratic in the fields and their canonical conjugate 
momenta) describing an infinite number of harmonic oscillators. Although these models can 
be solved, we will concentrate here on finite dimensional examples to avoid field theoretical 
complications, in particular, issues related to the presence of an infinite number of degrees 
of freedom and the coupling of the infinite different modes induced by the function /. 

We consider 

H = f(H ), H = ^ + tx 2 . (1.3) 

2m 2 

To simplify the calculations, we will assume that m = 1 and k = 1. We will also work with 
an arbitrary function / (subject to some mild smoothness conditions) until the very end of 
our discussion. At that point we will make some explicit calculations by using the functional 
form of the Einstein-Rosen Hamiltonian. We emphasize that similar arguments could be 
made for any system whose Hamiltonian is a function of a solvable one. 
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II. CLASSICAL TREATMENT 

We first discuss the classical solution for H = \{pq + xfy. The dynamics generated by 
H Q is given by the equations 

dx Q dp , . 

W =Po, ^ = -0. (2.1) 

Here we denote the time parameter as T because in the following we will compare two 
types of related dynamics where two different time parameters will be relevant. The general 
solution to these equations can be written as: 

x (T) = -^(ae~ iT + ae iT ), (2.2a) 
v2 

Po (T) = -|(ae-* T -ae ir ), (2.2b) 

where a and its complex conjugate, denoted as a, are fixed by the initial conditions (at 
T = 0) 

a = —jT- (2 ' 3) 

In view of this last expression it is useful to introduce a complex variable zq = xq + ipo to 
describe positions and momenta simultaneously. In particular Eq. ( 12.21) can be rewritten as 

z„(T) = z e- lT . (2.4) 

The trajectories in phase space, described now as the complex z-plane, are circumferences 
centered in the origin with radius \z \ = a/xq + pi = V^H . 

Consider next the solutions for H = f(Ho). To have well-defined equations of motion we 
require that the function / be differentiate. The equations of motion now read 

^ = {xJ(H )} = f(H )p, (2.5a) 
^ = {pJ(H )} = -f(H )x, (2.5b) 

where /' denotes the derivative of / with respect to its argument and { , } is the Poisson 
bracket. In principle, these coupled, non- linear, differential equations might seem difficult 
to solve, but because H is a constant of motion, 

we can simplify them by introducing a new time parameter 

T{t) = f'(H )t. (2.7) 
The reparametrization given by Eq. (12.71) allows us to transform Eq. (12.51) into the form of 
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Eq. (12.21) corresponding to the harmonic oscillator with the solution 



x 



(t) = x (T(t)) = [ae- i/,(Ho)t + ae l/ ' (Ho)i ] , (2.8a) 
v2 

p(t) = po(T(t)) = ^= [ae- lf ' {Ho)t - ae if ' {Ho)t ] . (2.8b) 
V2 

Note that although (x(t),p(t)) have the same physical meaning as (xo(T),p (T)), they are 
different functions - x(t) is the composition (in the mathematical sense) xo(T(t)) of xq(T) 
and T(t). We find in Eq. (12.81) an energy dependent definition of time that yields a different 
time evolution for each solution to the equations of motion. (H has a different value for 
each initial condition.) The orbits in phase space for H and f(H Q ) are the same taken as 
non-parametrized curves. Nevertheless, for H the curves are parametrized by T, whereas 
for f(Ho) they are parametrized by t. The solutions for Hq all have the same frequency 

z (T) = z e~ lT , (2.9) 

in contrast to those for f(H ) which have frequencies that depend on the initial conditions 
(through the value of H = \z \ 2 /2) 

z(t) = z e- itf ' {Ho) . (2.10) 

III. QUANTUM EVOLUTION 

The behavior of quantum systems is quite different from the classical one. We choose as 
a basis for the Hilbert space of the harmonic oscillator the states \n) which satisfy Ho\n) = 
h{n + l/2)\n). (In the following we choose units such that h = 1.) Every initial state can 
be expressed as 

oo 

|V) = 5>„|n>, (3-1) 

n=0 

and the evolution is given by: 

oo 

|* (T)> = e- li ^ T - To ^> = J2 c ne- l(T - ToKn+1/2) \n). (3.2) 

n=0 

Let us consider a Hamiltonian defined as H = f(H ). To define f(A) for a general self- 
adjoint operator A we must require that / satisfy the relevant conditions for the spectral 
theorems. 0] In our case any function defined on the spectrum of ifo would give rise to a 
well defined Hamiltonian, but because we want to discuss the semiclassical limit, we will 
require that / be differentiable. 

The eigenvectors \n) of the Hamiltonian Hq are also eigenvectors of /(-Ho) with eigenvalues 
f{n + 1/2). The evolution of a state |^) defined by f(H ) is given by 

oo 

|tf(t)) = e -if(Ho)(t-to)^ = Y^ Cn e- i{t - to)f{n+1/2) \n). (3.3) 

n=0 
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We see in Eq. (13.31) that the situation is not analogous to that found in the classical sys- 
tem. In the quantum mechanical case we cannot obtain from \^ (T)) by a simple 
reparametrization of time, even if we allow it to depend on the initial state vector 
because the relative phases between different energy eigenstates |n) change in time and pro- 
duce a non-trivial difference between the wave functions under the evolution defined by Hq 
and f(H ). 

IV. COHERENT STATES 

Once we know the exact classical evolution of any state, we can search for semiclassical 
states that evolve in the same way. In general, even for the harmonic oscillator, wave packets 
(more specifically their squared modulus) change shape as they evolve in time. 0, [1] However, 
there is a family of non- stationary coherent states whose wave function -0 (modulus squared) 
does not change its shape as time evolves. A plot of \ip\ 2 as a function of time shows that 
it rigidly moves back and forth as a particle subject to a restoring force proportional to the 
distance to a fixed point in space, that is, a classical harmonic oscillator with Hamiltonian 
H = \{p 2 + x 2 ) (see Fig. ED . 




FIG. 1: Evolution of the squared modulus of the wave function of a coherent state for the harmonic 
oscillator. 

These coherent states of the harmonic oscillator (and their free field counterparts) have 
a number of additional interesting properties including the following [2]: 

1. They are eigenstates of the annihilation operator -that can be written in terms of 
the position and momentum operators as a = -^(X + iP)- with complex eigenvalue z 
whose real and imaginary parts encode the initial position and momenta of the classical 
motion. In terms of z and its complex conjugate z we have (z\X\z) = ^(z + z) and 

{z\P\z) = 75(5 — z). If we start with the condition that \z) is an eigenstate of a, it is 
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straightforward to express \z) in terms of the energy eigenstates \n): 

00 n 

|,)=exp(-| 2 | 2 /2)J]^=|n). (4.1) 

n vn! 

2. The dispersion of the position and momentum operators in these states [(2;|X 2 |z) — 
(z|X|z) 2 ] = [(^l-P 2 !^) — (2|P|2;) 2 ] = | are constant and saturate the Heisenberg un- 
certainty inequalities (coherent states define minimal wave packets). It can be seen 
that coherent states are also minimal with respect to energy and momentum 0] in 
the sense that A Z H = \z\, the characteristic time of the system in the state \z) is 
t z = l/(2A g H ), and then t z A z H = 1/2. 

3. The time evolution of the state \zq) is given by 

e- lT "°\z ) = e- lT / 2 \e- iT z ). (4.2) 

This equation means that as time evolves, the unitary ray defined by a coherent state 
\z ) at T = (i.e. the set of vectors of the form e l9 \zo) with 9 6 R) remains coherent 
at any time T and is labeled by 

z(T) = e- iT z Q = x Q (T) + ip (T), (4.3) 

where the functions x$(T) and Po{T), given by Eq. f)2.2p . are the position and momen- 
tum of the classical harmonic oscillator as a function of time. 

4. The set of coherent states defines a linear, non-orthonormal, overcomplete basis of the 
Hilbert space for a harmonic oscillator. In particular, we can find a relation of the 
type 

- f d 2 ^)(^| = L (4.4) 

71" Jc 

As we can see, the coherent states for the harmonic oscillator satisfy a set of properties 
that allow us to consider them as semiclassical in the sense that their time evolution 
closely follows the classical one. They also satisfy interesting properties that render 
them an important tool in the study of oscillator systems or free field theories. 



V. EXAMPLE 

As an illustration of these methods, we will answer the question: Can we build appropriate 
coherent states for a one-particle system with a Hamiltonian of the form /(i?o) with Hq = 
|(p 2 +x 2 )? This case is interesting because if the answer were affirmative, it could be possible 
to extend the result for interesting field theories such as general relativity reductions of the 
Einstein-Rosen type. As we will see the answer to this question is in the negative. 

We show that it is not possible to build proper coherent states for f(H ) by proving that 
under time evolution the label z, which encodes the initial data, cannot evolve according 
to the classical dynamics dictated by f(H ). In terms of the initial data Zq = Xq + ipo 
(combined in the complex number z ), the classical evolution of the system is obtained from 
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Eq. (J2IIDD 



z(t) = z exp(-itf'(\z \ 2 /2)). 



(5.1) 



So we will require that the state, which we will also label \z) in analogy with the usual 
coherent states, evolve as 



exp(-itf(H ))\z) = exp(i<p(t,z)) \zexp(-itf(\z\ 2 /2))> . 



(5.2) 



Equation (15.21) is similar to Eq. (14. 2p . Note that we must work with unitary rays so we include 
an arbitrary phase exp(itp(t, z)). We now expand \z) in the orthonormal basis provided by 
the energy eigenf unctions |n) of the harmonic oscillator Hamiltonian 



(5.3) 



n=0 



where the coefficients ip n (z) = ip n (x + ip) = ip n (x,p) are taken as differentiable functions. 
Equation (15.21) becomes 



rp n (z) = exp[i(<p{t,z) + tf(E n ))]^ n (zexp(-itf'(\z\ 2 /2))). 
If we use the notation i(} n (x,p) for ip n (z) we can rewrite Eq. (15. 4p as 
ipn&p) = exp[i(<p(t,z)+tf(E n ))] 



(5.4) 



(5.5) 



x^nUcos(t/ / (|^|V2))+psin(tf(|z|V2)),pcos(t/ / (|z| 2 /2))-a;sin(t/ / (|z| 2 /2)) 



The left-hand side of Eq. (I5.5P does not depend on time, so the time derivative of the right- 
hand side must be zero. We evaluate this derivative at t — and obtain the consistency 
condition 



fW/2) 



p— — (x,p) — x— — (x,p) 



dx 



dp 



-i [<f (z) + f(E n ))ip n (x,p), 



(5.6) 



with (f (z) - Qt 
rewrite Eq. (I5.6P as 



d<p ( 0,z ' . By introducing polar coordinates x = rcosO and p = rsin^ we can 
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Equation (I5.7I) can be solved to give 

ip n {r,9) = c n (r 2 )exp 
where 

X(0,r) -- 



f(r 2 /2) 

f(E n )e + x{e,r) 



f'(r 2 /2) 



(5.7) 



(5.8) 



(5.9) 



and c n (r 2 ) are arbitrary functions of r 2 . It can be easily checked that for the usual harmonic 
oscillator, f(x) = x and E n = n + ~, the choice cp(t, z) = —t/2 gives ip n = c n (r 2 ) exp(in9). 
The latter can be written as ip n (z) = c n (|z| 2 ) exp(mArg n (2;)), where Arg n is a branch of the 
argument of z. With this choice exp(mArg n (2;)) is independent of the branch chosen for 
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the argument, and we can write ip n (z) = c n (\z\ 2 ) exp(inArg(z)) = c n (\z\ 2 )z n with Arg(z) = 
{arctan(jo/g) + 2kn,k G Z}. This result should be compared with the result ip n (z) = 
g-M l 2 z n /n\ corresponding to the harmonic oscillator coherent states. As we can see only 
part of the dependence on z is fixed by Eq. (I5.2p . but the result is compatible with ip n (z) = 
e H*l l 2 z n jn\. By using the other conditions the full dependence on z can be obtained. 

From Eq. (15.81) we observe that, in general, the result will depend on the branch chosen. 
This ambiguity is unacceptable so we conclude that it is usually impossible to have a family 
of coherent states satisfying a condition equivalent to Eq. (15. 2\\ for the evolution given by 
f(Ho). We consider an explicit example using the functional form given by the Hamiltonian 
of the Einstein- Rosen waves f(x) = 2(1 — e~ x ^ 2 ). The solution (15.81) for this choice is 

^ n (r, 9) = c n (r 2 ) exp (V 2/4 [2(l - e~^ n+ ^)9 + X{9, r)]\ . (5.10) 

We need to require that 

exp (V 2/4 [2(l - e-^ n+ ^)6 + X(9, r)]) (5.11) 

be independent of the branch chosen for the argument 9 -otherwise it is not single-valued. 
However, this requirement is impossible to satisfy because X(9,r) is independent of n. If 
we write 9 = 9 + 2kir, we obtain the condition 

4kne T ' 2/ \l - e-^ n+ ^) + e r2/4 X(9 + 2kn,r) =e r2/4 X(9,r) (mod27r). (5.12) 

If we consider Eq. (15.121) for two different numbers n and m, the difference between them 
gives 

4vrA;e r2/4 (e~^ n+ ^ - e^ (m+ ^ = (mod27r) (5.13) 
for all m and n which is impossible. 

VI. CONCLUSIONS 

We have described the classical and quantum dynamics of systems with Hamiltonians that 
are functions of other solvable Hamiltonians and compared the cases Ho and f(Ho). Classi- 
cally the states evolve in very similar ways and follow the same phase space orbits although 
with different time parametrizations. In contrast, their quantum evolution is qualitatively 
different due to the appearance of non-trivial relative phases. We discussed this issue by 
analyzing the existence of coherent states and their properties for functionally related Hamil- 
tonians. In particular, we gave a proof of the impossibility of constructing coherent states 
that satisfy the four conditions in Sec. IV for general Hamiltonians of the form H = f(Ho), 
with Ho corresponding to the harmonic oscillator. This case is especially significant because 
of the role played by harmonic oscillators in the description of free quantum field theories. 

By relaxing some of the conditions defining coherent states for the harmonic oscillator, 
we can conceivably find a set of suitable semiclassical states for the more complicated dy- 
namics given by f(H ). Our analysis does not exclude this possibility, but suggests that the 
definition that we must adopt will require major changes in the conditions that are satisfied 
by the familiar coherent states. 
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